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Abstract 

We study the motion of self deforming bodies with non zero angular momentum when 
the changing shape is known as a function of time. The conserved angular momentum with 
respect to the center of mass, when seen from a rotating frame, describes a curve on a sphere 
as it happens for the rigid body motion, though obeying a more complicated non-autonomous 
equation. We observe that if, after time AT, this curve is simple and closed, the deforming 
body 's orientation in space is fully characterized by an angle or phase 9m- We also give a 
reconstruction formula for this angle which generalizes R. Montgomery 's well known formula 
for the rigid body phase. Finally, we apply these techniques to obtain analytical results on 
the motion of deforming bodies in some concrete examples. 

MSC2000 Subject Classification Numbers: 53Z05 (Primary), 70F99 74A05 93B29 (secondary) 
JGP SC: Classical mechanics, real and complex differential geometry 

Keywords: Deformable bodies, reconstruction phases, time-dependent non-integrable classical 
systems 

Contents 

1 Introduction 2 

1.1 Background 2 

1.2 Main results |3 

2 Physical setting 4 

2.1 Deformable bodies 4 

2.2 Self deforming body hypothesis 5 

2.3 Equations of motion 6 

2.3.1 Gauge freedom |8 

* Corresponding author. E-mail: cabrera@mate.unlp.edu.ar 
Tel.:+54-221-4229850 int 105 Fax:+54-221-4245875. 



1 



Phases in the self deforming body motion 9 

3.1 Reconstruction 9 

3.1.1 Reconstruction of c(t) from c(t) in the bundle Qo — > Qq/SO(3) : S 

3.1.2 Reconstruction of R(t) from U(t) in the bundle 50(3) — > S 2 L : Jjj 

3.2 Generalized Montgomery formula 1C 

3.2.1 Preliminaries 11 

3.2.2 The formula \U 

Some Applications 

4.1 Solutions on the sphere 14 

4.1.1 Cases with I(t) = diag(h(t), I 2 (t), I 3 (i)) and L(d (t)) = for all i 15 

4.2 Examples LL7 



1 Introduction 

1.1 Background 

We are going to study the problem of describing the motion of a rotating body whose shape is 
changing with time in a known controlled fashion. A particular case of this problem is the one 
in which the body's shape is constant in time, i.e. a rigid body. 

As well known, a free rigid body rotates about its center of mass in a rather complicated 
way, depending on how its mass is distributed in space. This distribution is represented by the 
corresponding inertia tensor and the motion is such that the (spatial) angular momentum with 
respect to the center of mass is a conserved quantity. 

Analytically, the orientation of the body with respect to an inertial reference frame can be 
obtained by, first, solving Euler equations for the (body) angular momentum and, finally, re- 
constructing the desired curve in the space of rotations from the momentum one. A beautiful 
result by R. Montgomery ([3]) states that, when the momentum curve completes a period, the 
orientation of the body in space is the initial one up to a rotation in a certain angle about 
the (conserved) angular momentum direction. Moreover, he derived a reconstruction formula 
(see pQ) for that angle, usually called the rigid body phase, which involves a geometrical (an 
holonomy) and a dynamical (energy and period values) contribution. 

Now, when a rotating body is free, but not rigid because its shape changes with time in a 
prescribed fashion, the way in which the mass is distributed in space is thus also changing with 
time. How does such a body move? Or, since we know how its shape is changing: which is the 
rotation about the center of mass induced by this changing mass distribution? For self deforming 
bodies with zero angular momentum, this question was answered by Shapere and Wilczek in [8]. 
In that case, the induced reorientation has a pure geometric nature because it is described by 
a horizontal curve with respect to the mechanical connection in a S'0(3)-principal fiber bundle 
over shape space (see [8], [5] and references therein). 

Another related problem is that of finding the optimal sequence of deformations in order to 
induce a given reorientation of the deforming body. This is an optimal control problem which 
generalizes the well known falling cat problem (see [5]). The problem we want to analyze is, in a 
sense, the orthogonal to the above control problem: we know the sequence of deformations and 
we want to find the induced reorientation. 
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1.2 Main results 



In the present paper we shall focus on a case not covered in [Hj, i.e. the case in which a self 
deforming body rotates with non-zero (conserved) angular momentum. 

Our main result is an expression for an angle or phase that determines, at specific times, the 
exact orientation of a spinning self deforming body with non zero angular momentum, general- 
izing R. Montgomery's formula [I]. 

The examples that we shall be keeping in mind are the ones in which someone reaccommodates 
the furniture in a spacecraft or an antenna coming out from a satellite in orbit. 

Notice that in the above concrete examples, the body is acted by external forces (e.g. gravity). 
Nevertheless, also note that for small objects like satellites in orbit the angular momentum with 
respect to the center of mass is approximately conserved. Within this approximation, the full 
motion can be described by two sets of decoupled eqations: the ones for the center of mass (a 
central force problem) and the ones we shall give below for the rotation about the center of mass 
(a self deforming body problem). 

The total reorientation of a self deforming body has two contributions: the one induced from 
the change in its shape (of geometric nature [8]) and the one we shall study, that follows from 
having non vanishing angular momentum (of dynamical nature as for a rigid body). 

The class of deforming bodies we shall consider is the one that will be refered to as self 
deforming bodies. This bodies are defined by a pure kinematical constraint and a dynamical 
hypothesis described in section I2.2L In sec. 12.31 we shall derive the corresponding set of (second 
order) non autonomous equations of motion for the unknown rotation about the center of mass. 
These follow from the conservation of the angular momentum measured from a reference system 
having its origin at the center of mass and axes parallel to those of an inertial one for all time. 
We will refer to it as the spatial angular momentum. 

Also in 12.31 we shall observe that, as in the rigid body problem, the desired induced rotation 
can be reconstructed from a solution of the associated body angular momentum (first order, non 
autonomous) equations. This is the angular momentum as seen from a reference frame which 
is rotating with the deforming body (see [8]). At this point, we can re-state our main result: 
when, after some time, the body angular momentum solution returns to its initial value, the 
reconstructed rotation curve returns to its initial value up to a rotation about the (conserved) 
spatial angular momentum direction; moreover, in section [3] we show that the angle of this 
rotation or self deforming body phase can be expressed (mod. 2n) by the reconstruction formula 
(|13p involving a geometric and dynamic term. This result can be seen as a straightforward 
generalization of Montgomery's formula from the rigid body to the self deforming body motion. 

This formula relates the body's orientation with the (non conserved) energy integral and 
the geometry of the (non zero) body angular momentum solution curve. In the zero angular 
momentum case of [8], as the motion is of a pure geometrical nature, the above phase becomes 
trivial. 

As in the rigid body case, our formula can be applied when we have a geometric description of 
an underlying simple closed body angular momentum solution curve. Explicit time dependence 
of the equations implies that, in general, energy is not conserved during the motion of the 
body. Also, as the equations for the body angular momentum are non linear and have generic 
time-dependent coefficients, solutions are hard to describe in the general case. 

In view of this last observation, in sec. H]we complete this work by studying some particular 
classes of deformations. In each case, we shall be able to derive analytical results on the motion 
of the underlying deforming body by making simple dynamical estimates on the geometry of the 
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body angular momentum solutions and by thus applying the generalized Montgomery formula. 

Acknowledgements: A.C. would like to thank Dr. J. Solomin for stimulating discussions and 
suggestions. He would also like to thank CONICET- Argentina for financial support. 



2 Physical setting 
2.1 Deformable bodies 

Now, we review the setting presented in [5] (see also [8], [7]) for deformable bodies. 

Let us call Q the configurations space of a system of N— particles or an extended body 
from the reference system CM(t). Thus, Q = R 3Ar ~ 3 or Q = {smooth maps q : B C M 3 — ► 
M 3 s.t. f B dvol p{x) q(x) = 0} for B being a reference shape of the extended body. In both 
cases, the usual action of SO (3) on M 3 gives rise to an action of SO (3) on Q. This action turns 
out to be free on 

Qo = Q — Qid 

where Qid is the set of points in Q representing configurations in which all the particles or the 
entire body is contained in a straight line. Hence, 

Q A Q /SO(3) 

defines a principal fiber bundle, whose base B = Qq/SO(3) is usually called the shape 
space. 

In both particle system and extended body cases, the manifold Qo (and also Q) has a Rie- 
mannian structure induced by the usual scalar product of M 3 . So there is a natural principal 
connection on the bundle Qo — > Qo/SO{3) defined by choosing as the horizontal subspaces the 
orthogonal complement to the vertical subspaces with respect to this metric. This is usually 
called the mechanical connection on the bundle Qo — > Qo/ SO(3). 

Notation: From now on, 

• S will denote a given inertial reference frame, 

• CM(t) will denote the reference frame with origin at the center of mass of the 
body rcM(t) for each t and axes parallel to those of S, 

• CM(t) will denote any reference frame with origin at the center of mass of the 
body, with (possibly) rotating axes with respect to those of CM(t). 

Remark 2.1 (Reference systems) Note that a point qo G Qo over a shape Tt(qo) = bo G 
Qo/SO(3) gives the configuration of a body with shape represented by bo as seen from the 
reference system CM(t). Another point <?o s.t. ir(qo) = ^(qo) then represents the configuration, 
as seen from CM(t), of a body with the same shape but, now, rotated with respect to the one 
represented by qo- We can also interpret go as describing the same body but as seen from a ro- 
tated reference system CM(t). This last interpretation of the different points of a fiber 7r -1 (6o) 
is the one that we shall keep in mind for the rest of the paper. See also the discussion in ref. [8]. 
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2.2 Self deforming body hypothesis 



Let us call fi {t) the position of the i— th particle with respect to S(t) at time t. Then, for each 
time t, there exist a global rotation R(t) E SO(3) and a translation T{t) £ M 3 such that the 
position with respect to the inertial reference frame S is 



A self deforming body is defined to be a system of particles or an extended object satisfying: 

i) Kinematics: There exist a reference frame S(t), not necessarily inertial, from which we 
know fio(t) or, equivalently, a reference curve do(t) in Qq. Consequently, we also have a 
corresponding shape space curve c(t) = Tr(do(t)). 

ii) Dynamics: The constraint forces which act on the particles in order to give this prescribed 
motions fj (t) are internal forces satisfying the strong action-reaction principle. This 
means that all forces acting on the particle i are caused by other particles j 's and = 
—Fjf* with Ffj 1 * parallel to the vector = ri — rj. 

Condition (i) can be seen as a set of time dependent kinematical constraints generalizing the 
usual ones of rigidity: from S(t) we know how the body's shape is changing (see also [8]). 

Example 2.2 (Space-craft) For the system being a space-craft, S(t) could be chosen as a frame 
fixed to some part of the ship or an astronaut himself. 

Remark 2.3 (Mechanical forces) Notice that, although some forces do not satisfy the strong 
action-reaction principle (for instance, electro-magnetic forces), most of mechanical forces do. 

Remark 2.4 (Center of mass reference) We can always take S(t) = CM(t) (recall our notation 
12. If) having its origin at the center of mass at all time. See also the discussion at end of this 
section. 

Remark 2.5 (Non conservation of energy) Note that with these kind of time dependent con- 
straints, the energy is not conserved in general because the deformation is implemented by time 
dependent constraint forces. 

The self deforming body problem is to find a curve R(t) in 50(3) such that for 



in Qo the spatial angular momentum with respect to the cm. is conserved (see below). 
This can also be seen as a reconstruction problem (see pQ) for the rotation R(t) from the 
given c(t). 

We end this section with some remarks on the meaning and the measurement of do(t). First, 
we would like to stress that the reference curve do(t) is a natural physical input for the problem. 
To illustrate this fact, let us suppose that we want to describe the motion of a space-craft or 
satellite when someone is reordering the furniture inside of it, or when an antenna is coming out 
from this satellite. Before launching, in the lab., an engineer can attach the satellite to the floor 
and perform exactly the same deformation that will occur in space. The body does not rotate 
because it is attached, but the position of all its parts can be measured as a function of time t 



n = 



R(t)fio{t)+T{t). 



(1) 



c(t) = R(t) ■ do(t) 



(2) 
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from a lab. reference frame. Then, the position of the center of mass can be established for all t 
and, consequently, the position of every part of the body from a reference system CM(t) fixed 
to the center of mass can be known for each t. 

This provides us with a curve do(t) as desired: when the satellite is in orbit, the same defor- 
mation will occur yielding that do(t) projects onto the same curve in shape space as the physical 
curve c(t). Notice that, as the body can freely rotate about its center of mass, the position 
with respect to CM(t), represented by c(t), will differ, in general, by a rotation from the one 
described by do(t) for each t. This rotation is precisely the solution R(t) of the self deforming 
body problem. 

Example 2.6 (Rigid body) Note that the rigid body is a special case of the self deforming body: 
take fj (t) constant for all time. More generally, do(t) must be contained in the fiber over the 
point representing the constant shape of the rigid body for all t. 



2.3 Equations of motion 

The equations for R(t), according to our definition of the self deforming body, can be derived 
from the conservation of the angular momentum relative to the center of mass 

Lcm = 0. 

This means that the rotation must be such that, from a frame CM(t) this quantity is conserved 
even though things are moving internally in the system. 
Let us recall the well known quantities: for R(t) € 50(3) and d(t) = {ri(t)} £ Qq, 

• body angular velocity ~ R^R is defined, as usual, by x v = R~ 1 Rv for all 
v £ M 3 . We shall denote \& : (so(3), [,]) — ► (M 3 , x) the usual Lie algebra isomorphism 
(see for ex. [2]).; 

• (locked) Inertia tensor: I : Qq — » S^q 3 = {3x3 real symmetric positive definited 
matrices}, v ■ I({ri})w = Yl m i( v x r i) ' ( w x r i)'i 

i 

• Angular momentum (with respect to a rotated frame with origin at the center of mass): 
L : TQq -> M 3 , L({rj,rj}) = Y^ m i r i x h, satisfying 



L(£(R(t)d(t))) = R(t) I{d{t))"T + R(t) L{f t {dm 



(3) 



This gives the momentum map for the SO (3) action on TQq (see the details in [T], [5]); 
• Kinetic energy: T : TQq — > R, T({rj,rj}) = rf, for which 

i 

T(±(R(t)d(t))) = ■ I(d(t))"B {t) + LijtW))) ■ + T (j t ^)- W 

For the physical curve c(t) in Qq, the following quantity must be conserved: 

L CM = L{j t c{t)) = L(±(R(t)d(t))) = R(t) I(d (t))u;^ + R(t) L(±(d Q (t))) 
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Here I(do(t)) is interpreted as the inertia tensor measured from the reference frame S(t) = 
CM(t) and we shall call 

L (t) := L(—(d (t))) = ^2mir io (t) x r io (t) 

i 

the internal (or apparent [8]) angular momentum with respect to CM(t). 

The (time-dependent, second order) equations of motion for R(t) thus read 

j t L(R(t)d (t)) = (5) 

I(do(t))J B = I(d (t))LO B X uj b + L Q (t) x u B - ^-(I(d (t))) u B - 

dt dt 

when we express them in terms of the body angular velocity lob- 
The reconstruction equations for R(t), once we solved the previous one for lob, are 

R = R ujb (6) 

where Cjb = ^H^b) with ^ : (so(3), [, ]) — > (R 3 , x) the usual Lie algebra isomorphism (see 
[2]). The initial value R(ti) must be such that R(ti)do(t{) = c(t{) coincides with the initial 
value of the problem. 

Example 2.7 (Rigid body) For the rigid body, recall that do(t) must be contained on the fiber 
over a point in shape space. We can then choose the do(t) (equiv. fj (i)) to be constant for all 
t, so I(do(t)) = I is constant in time and L Q = 0. In this case, we recover Euler equations: 

Iwg = fog X 

as expected. 

a I r\ m r ttttI-K +riz^ t ln'in \r\r\r\-\T t -\ v «• \ W I aim a c / I = 



Also in analogy with the rigid body problem, as L(4c(i)) G M 3 is conserved during the time 



evolution, if we define 



Ii{t) = I(d,{t))o J f ) +L{j t d ) (7) 



we then have that L(-^c(t)) = R(t)U(t) and, hence, its M 3 — norm ||L(^c(i))|| = | ] IX 

is constant for all t. The quantity H(t) represents the angular momentum measured from the 
reference frame S(t) or body angular momentum. 

Remark 2.8 (Recovering the angular velocity) Since I(do(t)) is invertible for all i, we can 
recover at every time t the angular velocity u>b^ f rom H(t) £ M 3 : 
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(8) 



U = Ux(I-Hd (t))(U-L(ido(t)))) 



(9) 



Hih) = R- 1 (t l )L CM 

whose solutions lie entirely on the sphere 5£ C M 3 of radius \\L(^c(t)) 
Using ([5]), the reconstruction equations for R(t) become 



7T 



R = R V-^r^doWMt) - L(-d ))) 

dt 

or, equivalently, if we set R{t\) = id for simplicity 



(10) 



R(t) = Texp 




ds *-\i-\d (s))(n( s ) - LQdo(s)))) 



where T stands for the time ordered integral (see also [8]). 

Remark 2.9 (Non integrability) In general, as noted before, the explicit time dependence of 
the self deforming body tells us that energy is not conserved and consequently, we cannot 
reduce the dimension of the problem any further. 

2.3.1 Gauge freedom 

By definition, we are given a curve do(t) in the configuration space Qq, but we might want to work 
with another curve do(t) defining an equivalent self deforming body problem, i.e. 7r(do(i)) = 
w(do(t)) = c(t) G Qq/SO(3). This is equivalent to consider the self deforming body to be 

described from a new reference frame S(t) having the same origin and rotating, in a certain 
known way, with respect to the initial one S(t) from which the motion represented by do(t) was 
originally described. 

Remark 2.10 (Gauge transformations) This freedom in choosing the initial orientation curve 
do(t) can be seen as gauge freedom. Correspondingly, the change do(t) do(t) can be thought 
of as a gauge transformation. For more details on this analogy, we refer the interested reader 
to [5], and references therein. 

Among all possible lifts do(t) of c(t) we consider two: 

a) the horizontal lift with respect to the mechanical connection in the bundle Qo — ► 
Qo/SO(3). This is equivalent to the problem of finding a lift do(t) such that L(4do(t)) = 
Vt (see also remark \2. lip . 

b) a lift do(t) for which the inertia tensor I(da(t)) is diagonal for all t. This is equivalent 
to solve the problem of finding a lift of the base curve I(do(t)) along the map 



21 x 50(3) -► Slf 
(a,R) i — ► RaR- 

where 21 := {3 x 3 diagonal positive definited matrices}. 
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Remark 2.11 (Deformable bodies with zero angular momentum) In ref. [8], it is shown that, 
given a shape space curve, the motion of a self deforming body with zero angular momentum is 
described by the corresponding horizontal lift as in (a) above. These computations also arise in 
the falling cat problem (see [5] [7]) and other interesting problems (see references in [8]). 

Remark 2.12 (Simplifying the equations) Choosing a different do(t) changes the time-dependence 
of the coefficients of equation ([9]). Thus, an appropriate choice could turn this equation into 
a simpler equivalent one. For example, choosing the horizontal lift implies that the equation 
has vanishing L(do(t)) term because this is zero by construction. We also see that there is an 
obvious simplification when choosing the lift keeping the inertia tensor I(do(t)) diagonal. But 
observe that this two simplifications cannot always be carried out at the same time, since the 
horizontal lift does not necessarily diagonalize the inertia tensor in general. 



3 Phases in the self deforming body motion 
3.1 Reconstruction 

For completeness, we now describe two types of reconstruction phases QlJ) appearing in the 
configuration space during the motion of the body. In the rest of the paper, we shall focus only 
on the second (abelian) one. 

3.1.1 Reconstruction of c{t) from c(t) in the bundle Q$ — » Qq/SO(3) : 

Recall that, for each t, both do(t) and c(t) belong to the fiber over c(t) in shape space Qq/ SO(3) 
(see section I2.2[) . When the shape space curve is closed in [ii,i2]) we can then follow the 
standard procedure for reconstruction ([!]): choose do(t) to be the horizontal lift with 
respect to the mechanical connection having do(ii) = c(t\) as initial value. Then, do (£2) = 
Rg c (ti) with Rq being the holonomy of the base path c(t) measured from c(t\) with respect to 
this connection (see remark [2JJ] and sec. I2.ip . This is often called the (non abelian) geometric 
phase. Finally, under these assumptions, the reconstruction formula reads 

c(i 2 ) = R D (h) R G (h) c(ii) 

where Rd^) is usually called the (non abelian) dynamical phase. This dynamical phase can 
be obtained by solving eq. ([5]) with the initial value Ro(ti) = Id and with the above horizontal 
choice of do(t), i.e. with L(do) = 0. For details on general reconstruction see pQ. The interested 
reader can find details about this reconstruction for a deforming body motion with zero angular 
momentum in [8]. For a study of phases in the N = 3 body problem, we refer the interested 
reader to [6]. 

3.1.2 Reconstruction of R(t) from II(i) in the bundle SO (3) — > S 2 L : 

Recall that, in general, the unknown rotation R(t) in eq. ([2]) can be reconstructed via (|10p once 
we have solved the equation ([9]) on the sphere. An interesting special case is when this solution 
n(i) is closed in the interval [ti,t2] 5 that is when 

n(ti) = n(t 2 ). 
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In this case, there is a unique angle 9m naturally associated to this solution and to the initial 
condition R(t\) such that 

R(t 2 ) = exp(9 M pjj ) R(h), 

yielding 

c(t 2 ) = [exp(9 M -fa) R(h)] d (t 2 ) 

where L = ¥ -1 (.L) G so(3). We see that 9m defines an abelian reconstruction phase asso- 
ciated to the initial data R(t\) (coming from c(t\)). This phase appears when reconstructing 
R(t) from n(i) in a [/(l)-principal bundle 50(3) — ► 5^ that we shall describe in the next 
section. 

Remark 3.1 (Interpretation of 9m) Recall that R(t) takes the reference frame CM(t) to 
CM(t). This implies that, at time t 2 as above, the orientation of the body, as seen from CM(t), 
is precisely obtained by rotating the known configuration do(t 2 ) about the conserved angular 
momentum direction [Lqm) in the angle 9m- So this phase fully characterizes the position of 
the deforming body in space at specific times (i.e. t 2 ). 

In the rest of the paper, we shall focus on the latter reconstruction procedure. Note that, as the 
second phase is abelian, it is more likely to have simpler closed expressions for its reconstruction. 

Finally, we note that the most geometrically interesting situation is that in which both, the 
solution n(i) to Q and the shape space base curve c(t) in Qq/SO(3), are closed in the same 
interval [ti, t 2 ], i.e., 

n(ti) = n(ta) 

c(h) = c{t 2 ). 

When this conditions hold, there is a geometrically defined phase in the bundle Qq —>■ Qo/SO(3) 

c(t 2 ) = AR-cih) 

independent of the choice of do(t) (it only depends on the initial value c(t\)) and having the 
following expression: 

AR = exp(9 M |^) R^ARoR- 1 ^) , 

where ARq = Rq^^Rq 1 ^^ and R(t\) are fixed by the initial condition c(t\) and the angle 9m 
is again given by the generalized Montgomery formula presented in the next section. 

3.2 Generalized Montgomery formula 

In this subsection, we give a phase formula for the reconstruction of the rotation R(t) from a 
closed solution curve TI(i) of the equation Q. This formula generalizes the well known one given 
by R. Montgomery in [5j for the rigid body phase. For the proofs, we shall use some differential 
geometric results that we review below. 
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3.2.1 Preliminaries 



Recall the diagram (see [TJ) 

sol (3) A- T*50(3) ~ 50(3) x so* (3) sol (3) 

e<— too — »Adjj^ 

where: sol (3) denotes the Poisson manifold so*(3) with its (minus) standard Poisson bracket; 
7T and J are Poisson and anti- Poisson maps respectively and Ad* R := (Ad R -i)* denotes the (left) 
coadjoint action of 50(3) on so*(3). Recall (see e.g. [2]) that J is the momentum map associated 

to the left symplectic action of 50(3) on T*50(3). The trivialization T*50(3) ~ 50(3) x so* (3) 
by left translations is known as passing to body coordinates. 

If we fix an element L G so!_(3) ~ so(3) ~ M 3 (the isomorphisms are compatible with the 
corresponding Poisson brackets), then we have that 

*(Ad* R = R-*(0 

thus 

7r(J-\L))=S 2 L . 

The sphere 5£ of radius \\L\\ defines a symplectic leaf in sol (3) — so(3) ~ R 3 (see [2j). Moreover, 

in this case we have that 

J~ X {L) = {(R, II); R • n = L} ~ 50(3) Si 
(R,R~ l L) i — ► R~ X L 

is a [/(l)-principal fiber bundle (see [Tj). 

Now, consider the inclusion J _1 ( L ) ^ 50(3) x so*(3) ~ T*50(3) and the it(l) -valued 
1-form on J _1 (L) 

A := t^t-XQ L (11) 
WW 

where Q L is the canonical left invariant 1-form on T*50(3) in body coordinates. It can be seen 
that A gives a principal connection in the principal £7(1)— bundle J~ X (L) — ► 5£ ([TJ). This 
connection 1-form satisfies 

dA = -^-ri*u L 

WW 

where ui L = —dQ L denotes the canonical symplectic 2- form on T*50(3) in body coordinates. 
By the reduction theorem ([3], see also [2]), 

■* L * 
I UJ = IT U)^ 

with uifjL the reduced symplectic form on S\. Finally, if dS denotes the standard area 2-form on 
the sphere 5| C R 3 , then (see [2]) 
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3.2.2 The formula 



With these geometrical background, the following can be easily proved: 

Proposition 3.2 R(t) is a solution of the (second order) equation of motion Q iff (R(t) ,H(t)) € 
J _1 (L) C T*SO(3) is an integral curve of the time dependent vector field 

x(R,n,t) = (r *-\r\do(t))(n-L(do))), n x (/-Vo(t))(n-L(d )))). 

Remark 3.3 (Hamiltonization) This result can be viewed as a time dependent hamiltonization 
from TQq to T*SO(3) using the momentum map L and also a further reduction to J (L) of 
the problem equations of motion ([5|). See also similar comments about reduction for the 3-body 
problem phases in [6]. 

Thus, reconstructing R(t) from II(i) is the same as finding a curve (R(t),H(t)) G J~ 1 (L) 
in the U(l)— bundle J~ 1 (L) ~ SO(3) -—>■ S 2 L as above such that the projection to the base 
II(£) € is a solution of ([9j). Given n(i), we can apply the usual procedure of reconstruction 
([!]): choose Ro(t) € J~ l (L) in a natural geometric way as the horizontal lift of n(f ) € 
from the initial value -Ro(^i) = R(h) with respect to the connection A. Now, let 9(t) € be 
an angle to be determined by requiring the curve 

esp(0(t)pij) • (R (t),R^\t)L) = (exp(9(t)^-)R {) {t),R^ 1 (t)L) e J" 1 ^) - 50(3) 

to be the desired integral curve of X(R, II, t). In the above formula, L denotes ^~ 1 (L) E so(3). 
It follows that 9(t) must satisfy the following equation 

\\L\\9(t) = r\do(t))U(t)'U(t)-r\do{t))Lo(t)'U(t) (12) 

0(h) = o 

Now, note that if [ti, £2] C R is a closed interval, and II : [ti, t%] —>■ S\ is (any) continuous curve, 
then its image Im(Jl) is a compact, hence closed, subset of the sphere S\. So its complement 
Im{Ii) c is open and it exists a closed disc d C Im(H) c . We then have Jm(II) C dP and we thus 
showed 

Lemma 3.4 The image Im(JX) of a continuous map II : [^1^2] S\ is entirely contained in 
an open disc D C S\. 

We can now state our main result: 

Proposition 3.5 (generalized Montgomery formula): LetTl(t) be a solution o/Q satisfy- 
ing that II(ii) = II(i2) for some interval [ii, £2] that the image ofH : [ti, £2] —> is a simple 
closed curve (i.e. Imiji) homeomorphic to the circle S 1 ) then Rfo) = exp(9u jmr) R(h) an d 
the angle 9m is given ( mod 2tt ) by the formula 

(13) 

where D is a surface in S 1 ? bounded by the image ofTL. The — (resp. +) sign corresponds to the 
case in which the solid angle defined by D on the sphere, with its time-oriented boundary Ii(t), 
is a positive (resp. negative,) signed solid angle. 
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Remark 3.6 (Signed solid angles) As usual, we are considering a solid angle in the sphere to 
be positive or negative by applying the right hand rule to its oriented boundary (see [3]). Also 
notice that, mod. 2tt, the above formula keeps the same form (i.e., with the — sign) if we replace 

aT ^l^? by the corresponding signed solid angle. 

Remark 3.7 (Relation to the energy) The integrand in the right hand side of this formula can 
be expressed in terms of the total kinetic energy (see eq. (|1J)): 

(/- 1 (4)n(t) - j- 1 (d )^(4)) • n(i) = 

= 2T(j t (Rd )) - 2T(j t do) + r\do)L(d ) ■ L(d ) - r\d )L{d ) ■ U(t) 
Proof: By the above mentioned reconstruction procedure and since U(l) is abelian, 

R(t 2 ) = exp(e D ^) • exp(6 G jj^) • R(h) 

R(h), 

where 9n is the dynamical phase, solution of eq. (jl2l) and Qq the geometric phase, given by the 
holonomy of the base path H(t) with respect to the connection A and measured from R(tx). 
Thus the dynamical contribution Qd to 9m is precisely the second term in the r.h.s. of eq. (I13j ) . 

Let us then show that the remaining term coincides with the geometric contribution Oq. 

By the hypothesis and Lemma 13.41 Im(H) is entirely contained in a smooth disc D in S\. 
Being D contractible, the restricted principal [/(l)— bundle J~ 1 (L) \d — ► D is trivial and, then, 
we have a smooth section s : D — > J _1 (L). Once we have chosen the disk D containing the 
curve Im(U), the existence of a surface D C S\ whose boundary is Il(t) is obvious since D is 
diffeomorphic to an open disk in M? and Im(U) is homeomorphic to S . Thus, mod 2ir, we can 
write (see [I]) 

G = - I ! s*(dA) 




D 



ds- area{D) 



||i/|| J JD \\L\\ 

when the solid angle defined by D is positively oriented with respect to the (time oriented) 
boundary curve 11(f). The last two equalities follow from the results reviewed in the previous 
section. Formula fl 13 j) is then completed. □ 



Example 3.8 (Rigid body) For the rigid body, the kinetic energy T is conserved and, as we 
observed previously, do(t) can be taken as a point for all t. So L(J^do) = and the inertia tensor 
I (do) = I is constant. In this case, the periodic solutions of Euler equations bound disks on 
the sphere, thus D defines the usual signed solid angle and the above formula becomes the well 
known reconstruction formula derived by R. Montgomery ([!]). 
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4 Some Applications 



4.1 Solutions on the sphere 

We shall now describe some tools which can be used to study the geometry of solutions of eq. ([9]) 
on the sphere S\. Focusing on some particular cases we will be able to use this characterization 
of the solutions to yield analytical results on the motion of self deforming bodies by calculating 
the associated generalized Montgomery phase 9m- 

• Reconstruction of R{t): When the solution IT(t) for some time interval [£a,£,b] is an 
open path, we noted before that the rotation R(t) can be expressed as exp(9(t)jj^ ) • 

(Ro^),^ 1 ^)^), with (Ro^t),^ 1 ^)^ the horizontal lift of the base path U(t) with re- 
spect to the connection (JXTJ) and 9(t) a solution of eq. (fl2j) . When the solution IT(t) is 
a closed simple curve for a time interval [t^s]; we have a well defined phase 9 m given 
by formula (|13p . So, given a solution IT(t) in [ii.^L we can find a total phase by adding 
phases corresponding to sub-time intervals fa ti+i] for which the solution is a simple open 
arc in S? or a closed simple curve in S%. In the first case, we have a phase defined by 

R(U+i) = exp(9(t i+1 )j—)Par(R(ti)) 

with Par : 7r _1 (n(tj)) — ► 7r~ 1 (n(tj + i)) the parallel transport (see PQ) in the U(l)— principal 
bundle J _1 (L) — S\ of the initial condition R(t{), and 9(t) the solution of (|12|) with 
9(ti) = 0. In the second case, fixing the initial value R(ti), the phase is defined by 
R(ti + i) = exp(0M L) R(U) with 9m given by formula fjT3|) . 

• The Energy: As we noted before, in general, the energy is not a conserved quantity 
during the self deforming body motion. Nevertheless, if we know the evolution of the 
kinetic energy T(f t (Rd )) with time, we will be able to determine a specific subset of S£ 
in which the corresponding solution IT(t) lies. This fact can be shown as follows: let us 
define for each time t 

E t : S 2 L ^R 

: u^^u-r\d (t))u. 

Note that 

E t (U(t)) = T(±(Rd )(t)) - T(|db(t)) + \L(do(t)) ■ r\d (t))L(do(t)) 

for Il(i) a solution of ([9]). Hence, as do(t) is given, E t {U.{t)) is uniquely determined by the 
kinetic energy T(4?(Rdo)(t)). In this case, the corresponding solution H(t) on the sphere 
at time t must lie in the set 

E t -\k(t))ns 2 L 

where 

k(t) = T(j t (Rd )(t)) - T(j t d (t)) + h(d (t)) ■ rHdoimidoit)). 

The level sets E^~ 1 (k(t)) are (generally non centered and rotated) ellipsoids for each k ^ 
and each t. Also notice that, for a fixed time tj, the intersection E^. 1 (k(ti)) n S\ gives the 
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set where the body angular momentum of a rigid body with constant inertia tensor equal to 
I(do(to)) an d energy k{t\) would lie. Finally, the equation for the evolution of Et(H(t)) is 

jEmt)) = P(t) x i-\d (t))u(t)] ■ r^Mmm)) + in(t) • j t [i-HMt))Mt) 

which is coupled to the equation ([9]) for II(i). 

The arc-length: for a given closed time interval [ti, ^2] we are going to find a bound for 
the length of II([£i , ^2])- To that end, we note that 



d „, . 
dt^ 



uxir^doimu-Lidoit)))) 
< ||n||(||/- 1 (do(t))n|| + |/- 1 (do(t))L(do(/))) ) 

and, if p' 1 {d (t))v\\ < a _1 (t) \\v\\ for all v E M 3 , t € [h,t 2 ], then 

< ||n|| 2 a" 1 ^) + 1 1 XX 1 1 a" 1 ^) L(do(*))) 



Since llnl 



|L|| = I is constant, we then have that 



lenght(n([ti,t 2 })) < I / a~ 1 (t)(Z+ L(d (t)) )dt 



When a _1 (t) is a very small function (compared to y), we can deduce that II([£i, t 2 \) 

is contained in a small patch in ££. 

For general time dependent parameters I~ 1 (do(t)) and L(do(t)) we cannot give a characteri- 
zation of the solution II (t) of eq. ([9|). So we shall focus on some specific cases to illustrate how 
to handle concrete problems. 

4.1.1 Cases with I(t) = diag(Ii(t),h(t),I 3 (t)) and L(d (t)) = for all t. 

Let us denote by (1,2,3) the cartesian axes of so*(3) ~ M 3 and, hence, II = (LTi , II2 , II3) £ R 3 . 
In this cases, the intersections of each axis with the sphere S\ give a constant solution of (|9|), 
because at that points II is parallel to J _1 (i)II and the r.h.s. of eq. ([9|) vanishes. 
The equation for the evolution of the energy becomes 

±E t m)) = \m-j t [i-\d (t))}m) 



and the arc-length is bounded by 

lenght(U([t 1 ,t 2 ))) < I 2 a -1 (t)cft. 

Jti 

Now, we shall analyze further the special case 



h(t) < hit) < I 3 (t) 



15 



for all t E [iijfa]- Notice that this is the case (up to renumbering the I{s) for small enough time 
intervals [£i,£2]- Under this conditions, the axes of the ellipsoids E^ l (k(t)) coincide with the 
cartesian axes in M 3 and the arc-length is thus bounded by 



lenght(U([ti,t 2 ))) <l 2 I^ifdt 

Ju 



Fixing the time t, we have that through each point of S\ passes a solution of Euler equa- 
tions (rigid body) with inertia tensor diag(I\(t) , 12(f) ; I$(f)) ■ F° r each time we then have the 
corresponding homoclinic solutions (see for ex. given by the intersection of S\ with the 
ellipsoid of energy k(t) = j^y- 

Given a solution II(i) = (Ili(i), H 2 (t), II3 (t)) of ([9]) for the interval [£1, ^2] with initial value 
II(tx), the function f(t) = Et(IL(t)) reaches a maximum and a minimum on [ti, £2]; denoted E max 
and E m i n respectively. The same happens with the value of the principal moments of inertia 
Ii(t). The solution II(t) is then contained in a connected "crown like" region R which is the 
connected component of U te [ tl ^ 2 ]E^~ ([E m i n , E max ]) n S\ which contains the initial value H(ti). 

We can now show the following results on the qualitative behavior of H(t): 



1. if 



h 



then R is contained in either the semi-space IT > or in IT < 0. In this 



case, II(t) evolves in S\ describing a trajectory that orbits surrounding the cartesian axis 
1. More precisely, if the initial point lies in the component with (say) IT > 0, then the 
solution will lie in this component for all t in [ti, £2]- So if we consider spherical coordinates 

111 = I cos6 

112 = I sinQ coscp 

113 = I sinO sirup 

with 9 G [0,7r], (p € [0, 2ir], it follows that 6(f) < \ for all t in [ii,^]- From the equation 
([9]) we can deduce that 

j t <p = cos6 [I^(t) - I^(t) + (IzHt) - I^{t)) sin\\ (14) 

thus, being I\(f) < I 2 (f) < 13(f), then ip(t) is a monotonous decreasing function of time, 
showing that the solution II(£) tends to describe revolutions about the 1 axis. 



2. if 



then R is contained in either the semi-space II3 > or in II3 < 0. In 



this case, TI(£) evolves in S\ describing a trajectory that orbits surrounding the cartesian 
axis 3, as in the previous case. 

3. In other cases, the solution can pass from orbiting one axis to orbit another one. 
To show this, let us suppose that E^ili^i)) > — and that I\ is constant. Then 

lenght(U([t 1 ,t 2 })) < l 2 Ii l (t 2 - h) and so we can choose I\ such that n([*i , £2]) is con- 
tained in some small patch in S\. In the case that I 2 is also constant in time and I$(f) 
grows (note that the order is maintained in time), then 4fE t (U(t)) = ^-3(^)^3 (t) < 0. 
So the energy decreases as fast as we want if we make I$(t) grow sufficiently fast. Note 
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that n|(t) is bounded from below because II([ii,i2]) is in a small patch. In this situation, 
Et 2 (TL(t2)) can be made smaller than j-* — , so the solution is able to pass from the regime 
(1) to the regime (2) described above when the energy Et "crosses" the homoclinic energy 
boundary ^ — . 

-* 2min 

Remark 4.1 (Return time) In either of the previous cases (1) or (2) we can give a lower bound 
for the (shortest) return time T = t2~ t\ s.t. II(f i) = II^)- If we suppose that the solution 
starting at II(ii) satisfies the conditions of (1) above and that it returns to this value for the 
first time at t2, then 

T = , _, > 2^1^) x (1,0,0)|| 
2 1_ PIT 1 

±max 

The case corresponding to (2) is analogue. 

Now, suppose that we are in the case considered in (1) ((2) is analogous) above and that, 
in some interval C [ti,t2], the solution H(t) describes a closed simple curve in S\. 

Then, we can apply formula (|13f) to find the corresponding phase. Taking into account the time 
orientation of the closed solution II(i) (fixed by (|14p ). if are ^(- p ) < 2ir (resp. > 2ir) we must then 
take the + (resp. — ) sign in (fT3|) and we have that 

area(D) 2 area(D) 2 

± p 1" J Emin [U+l — H) < "M < ± p ^ J ^max {U+l ~ Uj- (15) 

Remark 4.2 (Bounding 9m mod. 2ir) Note that, since 9m is defined mod. 2tt, the above 
bounds yield nontrivial information when | (E max — E m i n ) (U+i — t{) < 2ir. 

4.2 Examples 

We now apply the previous techniques to obtain estimates for the motion of simple classes of 
deforming bodies. 

Example 4.3 (Vibrational deformation) In this case, we suppose that the body is globally 
shrinking or expanding, that is, the position of a particle from the given reference frame S is 



r io (t) = a(t)r, 



to 



where rj is a constant vector and a(t) is a never vanishing positive scale factor. This means 
that we can choose the curve do(t) in Q such that 

I(d (t))=a 2 (t)I 

with Ig the constant inertia tensor corresponding to the constant configuration {ri }. By a 
constant rotation, we can choose the reference system S (equivalently, another curve do(t)) from 
which Iq is diagonal. Then, eq. Q on the sphere becomes 



^n = a- 2 (t)(nx/- 1 n). 



Given an initial value II(ti = 0), eq. © can be exactly solved yielding 



u(t) = u RB ([ t 

Jo 



a (s) ds), 
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/(n) = -n-/ - 1 n 



with Urb denoting the rigid body solution of Euler equations IT = IT x I Q 1 H with initial value 
IT(ti = 0). Now, the function 

1 

2 1 

is constant along the solutions. Note that II(i) describes a closed simple curve on the sphere for 
t € [0, T] when J Q T a~ 2 (t) dt = Trb equals the period of the rigid body solution Urb- In that 
case, the corresponding phase is 



9m = -Arb + ttIt, /(IT) f a- 2 (t) dt 

wm Jo 



= -Arb + ijljr /(n)TRB 

where Arb is the (signed) solid angle enclosed by the rigid body periodic solution Urb with 
energy /(IT). Notice that this phase coincides with the rigid body phase for Urb ([4J). The motion 
of this kind of vibrating bodies is similar to rigid body motion up to a time reparameterization 
which is induced by the expansion/contraction. 

Example 4.4 (Expansion/ Contraction of an axially symmetric body): Let us consider the case 
of an axially symmetric body which expands in the direction of its symmetry axis, i.e., the case 
in which there exists a curve do(t) such that 

I{d {t))=diag{I 1 (t),l2,h), 
with I2 = Is- As in the previous case, eq. Q can be exactly solved: 

n(t) = n RB ( f ( J r'(*0- J 3~') ds) 

with Urb the rigid body solution to Euler equations IT = IT x / _1 ((io(ii = 0))IT and initial value 
H(t\ = 0). The function /(IT) = ^IT • /~ 1 (0)II is again constant along the solution IT(£), which 

is closed simple curve for t G [0, T] when J*„ T ^-1^! ^-1^ ds equals the rigid body period Trb 

^1 [0)— 1 2 

corresponding to Urb- In that case, the associated phase is 



Arb + -Jrrr f U(t) ■ I- l (do(t)) U(t)dt, 
\\ L \\ Jo 



where Arb is the (signed) solid angle enclosed by the rigid body periodic solution Urb with 
constant energy /(IT). Notice that, in general, this phase is different from the rigid body phase 
associated to iT^- 

Example 4.5 (An antenna coming out from a satellite along a principal axis): We now consider 
the cases in which 

HMt)) = diag(h(t),h ,h ) 

or 

I((k(t)) = diag(I lo ,I 2o ,h{t)) 

with Ii(t) (or Is(t)) an increasing function of time. These cases give simplified models for the 
situation in which an antenna comes out from an orbiting satellite along one of the principal 
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axes of inertia 1 or 3. Note that the satellite is free to rotate around its center of mass and so 
its motion can be described by eq. dSJ). Suppose that initially 1% < I2 < h- Then, in the first 
case, as Ii(t) grows this relation might stop holding after some time, so the solution could pass 
from orbiting one axis to orbit another one. Consequently, We have no control on this kind of 
solution. More precisely, as 

is negative implying that the energy decreases and the solution can pass from the case (1) to (2) 
of the previous section, describing an open curve on the sphere which we cannot characterize in 
general. In turn, in the second case the ordering prevails and 

is also negative. Since the energy decreases, if the initial value II(ti) corresponds to case (2) of the 
previous section, the solution also evolves according to (2) and we have a good characterization 
of its behavior. In particular, if n([ij, ij+i]) is a closed simple curve, we then know that the 
corresponding reconstructed rotation R(ti+\) is exp(0M jj^) R(ti) where from (fToll . 

area(D) 2 area(D) 2 

— T ill ^ TTFTT ^ min V'i+l ~ ^i) S "hi S ± — 2 ^ 17771 ^initial \ti+l ~ ti), 

\\L\\ W-L'W \\L\\ 



with Ei n itiai the initial (hence the maximum) value of the energy E t in [ti, ti + ±] and E m i n = E ti+1 . 
We thus note that we can have a better description of the motion of the satellite when the antenna 
comes out along the largest principal axis of inertia. 

Remark 4.6 (Slow deformations) Intuitively, when the antenna comes out very slowly, the 
motion of the satellite will be close to a rigid body motion. This is reflected in the fact that 
when ^[/jT 1 ^)] is very small (compared to trien E min ~ E in i tia i and "^jf - ~ A-RB- 

So the phase 6m is approximately the same as the rigid body phase associated to a rigid body 
with I = /(do(tj)) and initial II(tj). 

Remark 4.7 (Small bodies in the gravitational field) For a body, v.g. a satellite orbiting the 
earth, which is small with respect to the interaction distance with another body (eg: the earth), 
it is a very good approximation to suppose that the gravitational force acting on a particle of 
mass m.j of this body is 

-G M . 
F ' = m (r CM - Pf ^ " } 
where rcu denotes the position of the center of mass of the body. P and M denote the position 
of the center of mass and the total mass of the second body (eg: the earth), respectively. We thus 
see that the equations of motion for the position of the center of mass (a central force problem) 
are totally decoupled from the equations of motion giving the rotation about the center of mass 
(a self deforming body problem as in ex. I4.5[) . 
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